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Abstract 

The dyadic Green functions for a uniaxial dielectric-magnetic medium, together with a re- 
versible field transformation, were implemented to derive a formulation of the Huygens princi- 
ple appropriate to a uniaxial dielectric-magnetic medium with gyrotropic-like magnetoelectric 
properties. 
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1 Introduction 

Mathematically, the Huygens principle is represented by a formula in which the frequency-domain 
electromagnetic fields in a closed, source-free region are expressed in terms of the tangential com- 
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ponents of the electromagnetic fields on the surface bounding the region. A convenient means of 
expressing this formula is provided by the dyadic Green functions appropriate to the medium filling 
the region (Felsen & Marcuvitz, 1994). Accordingly, the Huygens principle plays a central role in 
electromagnetic scattering theory. It may be applied, for example, to the analysis of diffraction 
from an aperture, wherein the aperture is formally represented as an equivalent source (Chen, 
1983); likewise, it can be used to formulate the Ewald-Oseen extinction theorem and the T-matrix 
method (Lakhtakia, 1994). 

While the Huygens principle for isotropic dielectric-magnetic mediums is a staple of standard 
textbooks, exact formulations for more complex mediums are generally unavailable. This general 
unavailability mirrors the scarcity of exact formulations of dyadic Green functions for complex 
mediums. Notable exceptions to this include the Huygens principle for isotropic chiral mediums 
(Lakhtakia, 1992) and uniaxial dielectric mediums (Bergstein & Zachos, 1966; Lakhtakia, Varadan, 
& Varadan, 1989). A formulation of the Huygens principle for a generally anisotropic dielectric 
medium was established by Ogg (1971), but it does not yield explicit expressions as it employs 
a spectral representation of the dyadic Green functions. Tai (1994) developed a formulation of 
the Huygens principle for a simply-moving isotropic dielectric medium, based on the Minkowski 
constitutive relations. However, these constitutive relations describe only nondispersive mediums 
and are less general than the ones implemented in this communication. We present here a for- 
mulation of the Huygens principle for a general uniaxial dielectric-magnetic medium (Weiglhofer, 
1990) with gyrotropic magnetoelectric properties (Lakhtakia & Weiglhofer, 1990). A harmonic 
time-dependence of exp {—iwt) is implicit in the following sections. 

2 Preliminaries 

Let us consider a homogeneous bianisotropic medium characterized by the frequency-domain con- 
stitutive relations 



D{r) = I • E{r) + T x H{r) 
B{r) = -r X E{r) + M • H{r) 




(1) 



where the permittivity and permeability dyadics are 
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(3) 



with 7 being the 3x3 identity dyadic^ . Hence, the medium exhibits uniaxial dielectric and uniaxial 
magnetic properties, with the unit vector u indicating the direction of the distinguished axis, 
and the gyrotropic-like magnetoelectric properties mediated by the vector F. The form of the 
constitutive relations (1) is sufficiently general that a wide range of natural and artificial uniaxial 
mediums are described by them (Mackay & Lakhtakia, 2008). The incorporation of gyrotropic- 
like magnetoelectric properties extends the range of mediums described by eqs. (1) to include 
exotic metamaterials as well as certain simply moving mediums (Lakhtakia & Mackay, 2006a) and 
mediums inspired by Tamm's noncovariant formulation for vacuum in general relativity (Lakhtakia 
& Mackay, 2006b). 

The frequency-domain Maxwell curl postulates can be set down as 

V X E{r) + icoT x E{r) - iuifi • H{r) = -J_m{r) 

V X H{r) + iujT X H{r) + iujg • Mx) = leiL) 
wherein J g and are arbitrary electric and magnetic source current density phasors, respectively. 
In order to simpliiy the analysis, we use the transformations (Lakhtakia & Weiglhofer, 1997) 

e(r) = E_(r_) exp (iwF • r) 
h{r) = H_{r) exp (iujT • r) 
iJ(L) = Jeir) exp (iujr • r) 
l^iz) = Imir.) exp {iojT • r) ^ 
of the field and source phasors. Thereby, eqs. (3) transform to 

V X e(r) - iujg • h{r) = -J „(z:) 

V X h{r) + ioje • e(r) = j ^(r) 
which yield the Helmholtzian equations 

(V X X) • • (V X X) - • e(r) = iuije{r) - (V x X) • /x"^ • j (r 



(4) 



(5) 



(V x|) 'g-^ • (V xX) -w^^ 



3 Dyadic Green functions 



h{r) = iojjmir) + (V x X) 'g ^ • i^{r) 



(6) 



Let us now reproduce the dyadic Green functions for a uniaxial dielectric magnetic medium, as 
derived by Weiglhofer (1990), for convenience as they are essential for the analysis in the next 
section. 



^Here and hereafter, vectors are underlined whereas dyadics are double underlined. 
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The linearity of eqs. (6) ensures that their solution may be expressed as 

h{r)=[[[\g {r,r') • j (r') + g {r,r') • j Jr')] d\' 



(7) 



(8) 



where the integrations range over the region where j ^ ^ are not null valued. Herein, the 3x3 
dyadic Green functions g (a/3 = ee, em, me, mm) satisfy the equations 

= a/3 



• .9 (ziiZl') + V X g {r,r') — £S{r — r') 

— =ee = me — 

V X g {r,r') — iojfj. • g {r,r!) = 

= ee = =me ~ 

V X £f (r , r') + iuje • g (r , r') = 

= mm — =em — 

—iuiiJ.'g {r,r')+\^xg ir,r') = —]^S{r~r') 

= =mm =em ~ 

which decouple to provide the system of second order differential equations 
(V X X) (V X X) -w^i .g {r,r') =iu!lS{r-r') 

(V X X) • • (V X X) - cj2| .g {r,r') = - (\7 X I) ' fjr^ 6{r- r') 
(Vxl) 'I"! • (Vxl)-wV] 'g (r,r') = (V xl) •|-M(r-r') 



(V xl) 'I"! • (Vxl) 



(9) 



g iZiZ') = iw/^(r — r') 

= mm — 

with S representing the Dirac delta function. These are the dyadic Green functions of a uniaxial 
dielectric-magnetic medium which are known to be (Weiglhofer, 1990) 

VV 



9 {r,r!)=iLLiiJL 



_T(r-r')+ e„e-i 



ge{r-r') 



9 (r,/) = -MM"' • (Vx7) • [-e„5e(r:-r')e-i +T(r:-r')] 

= me = ^ — ' '- — — 

VV 



9 (H) e') = e 



T(r-r')+ M«M"' 



gt,{r-r') 



(10) 



The scalar Green functions g^^ in eqs. (10) are defined by 



exp 



9r,{R) 



47r(^7?„E«^-i 'R) 



1/2 



(»? = e,M), 



(11) 
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while the dyadic function T is specified as 

{Rxu){Rxu) 



m) = 



{Rxuy 



e jj, 



+ 



ge{R) {euR ' • E)'^' - {f^uR ' • r) 
iu) (e/x)^^^ {R X u)^ 



1/2 



I — uu ■ 



2{Rxu){Rxu) 



(12) 



{Rxuf 

On considering the dyadic transposes, which we denote by the superscript 'T', eqs. (10) reveal the 
symmetries 



9 iZiZ') = 9 {r!,r) 

= ee =ee 

= em =me 

= me =em 



(13) 



while interchanging r and r' yields 



5 {L,r')=g {r',r) 

= ee = ee 

.9 fcz:') = -5' (r'iZ:) 

= em = em 

5 fciO = (i:',2:) 
9 ir,r')=g {r',r) 



(14) 



4 Huygens principle 

The dyadic Green functions (10) may be exploited to derive mathematical statements of the Huy- 
gens principle pertaining to radiation and scattering in a uniaxial dielectric-magnetic medium with 
gyrotropic-like magnetoelectric properties. Consider a source-free region Ve, which is occupied by 

a medium described by the constitutive relations (1) and enclosed by the finite surfaces S and 6*005 
as schematically illustrated in Fig. 1. The source phasors J g „ (and j are nonzero only inside 
S. The unit vector n on 5 U Soo is directed into Ve- 
We introduce the dyadic function 

A{L,r') = -e{r) X J^lf^ ' S/xg (r,r') | - • [V x e(r)]| x 5 (r,r'); (15) 
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and note that 



V • A{r,r') = e(r) • (v x |/£"^ • V x g {r,r') }) " x {/£"^ * [V x e(r)]|) • g {r,r'), 

(16) 

by virtue of the identity 



V' (px£) = (Vxp) 'g-p' (Vx£), 



(17) 



and the fact that the dyadic /i ^ is symmetric. The dyadic divergence theorem (Tai, 1994; van 
Bladel, 1985) 



V - A{r,r') (fr 



Ve 



SUSo. 



(18) 



then yields 



(r) • X l^i^i • V X £ (r,r') }) " (v x j^'^ • [V x e(r)]|) • g {r,r') (fr 

= j ^ [ [^(H) X e(r:)] • {^-1 • [V xg^^(r,r')]} 

+ (^(e) X 1^"^ • [V X e(r)]|j (r,r') dV, (19) 

after using the identity s • x = (s x p) • Q. Now eq. (6)i and (9)i can be employed, 
together with the symmetry of the constitutive dyadic e and the fact that there are no sources in 
Ve, to evaluate the integral on the left side of eq. (19). Thus, we find 



ioje{r') 



[n(r) X e(r)] • \ /i ^ • V x g (r,r') > 

SUS^ L L= L =ee J J 



+ (^^(z:) X • [V X e(r)]|^ • g {r,r') (fr, r' G Vg. 



(20) 



Next, let 5*00 be sufhciontly distant from S that the integral herein on the surface Soo is eliminated 
by virtue of satisfaction of appropriate radiation conditions by the dyadic Green fimctions (Fclscn 
& Marcuvitz, 1994). Equations (5)i and (8)2 allow the curl terms in the integrand of eq. (20) to 
be replaced by explicit terms; hence, 

e{r')= [ f \\n{r) X e{r)] ' g {r,r') + \n{r)xh{r)]'g (r,r')| rfV, r' eVe- (21) 

J Js =me =ee ) 

On interchanging r and r', and exploiting the symmetries of g , g and g listed in eqs. (13) 

= ee = em =me 

and (14), it emerges that 

e(r) = / [ \ g {r,r!) ' [n(r') x h{r')] - g {r,r') • [n(r') x e(r')]| dV, L G (22) 

J J g K=ee = em J 



And, lastly, after inverting the transformation of presented in eq. (4)i, the Huygens principle 
for the electric field phasor emerges as 

Ml) = f fla ' [n{r') X Hir^)]- g {r,r!) ' [nir!) x E{r')]] 

J J V'=-ee — em J 

X exp [-iuiT ' (r - r')] d^r', r G (23) 

To find the corresponding expression for the magnetic field phasor, we take the curl of both 
sides of eq. (22) and utilize eqs. (8)2,3 to get 

V X e(r) = / / ( it^M * .9 {LiI!) ' [^(n') x 
J Js ^ = 

+ X |e-^ ' IVxg {r,r')] ] • \n{r') x e{r')]) d'^r', r G V^. (24) 

lU) L~ L =mm J J / 

Then, exploitation of eqs. (5)i and (9)4 to eliminate the curl terms delivers 

h{r) = (jig (r,r') • [n(r')xft(r')]-5 (r,r') • [n(r') x e(r:')]| dV, 

J Jq L =me =mm ) 

r G r ^ 5, (25) 

wherein the observation point r is explicitly excluded from the surface of the source region. The 
Huygens principle for the magnetic field phasor thus emerges as 

Eir) = [jig (r,r')- [n(r')x^(r:')]-5 (l, r') • [n(r') x £(r:')]| 

X exp [-iojT • (r - r')] d'^r', rGVe,r^S, (26) 

upon inverting the transformation of presented in eq. (4)2. 

5 Closing remarks 

The main results of this communication arc eqs. (23) and (26), which represent formulations of the 
Huygens principle appropriate to the bianisotropic medium described by the constitutive relations 
(1). This is the most general formulation of the Huygens principle presently available, as far as we 
are aware. We note that the formulation of the Huygens principle for a uniaxial dielectric medium 
derived by Lakhtakia, Varadan, & Varadan (1989) follows immediately from eqs. (23) and (26) 
upon setting F = and = /ioX , where /lo is the permeability of free space. 
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